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The most recent experimental results on the cosmic ray composition and energy spectrum above 1 TeV are 
reviewed and discussed. All data agree on the presence of the so-called "knee" at an energy 
the knee is seen in all the components of the Extensive Air Showers. These results support the hypothesis of an 
astrophysical origin of the knee, while no new features in the hadronic interactions at high energies are envisaged. 
The cosmic ray composition below and above the knee region is still an open question. According to most 
experiments, the knee seems to be due to the light component of the primary beam, with a composition getting 
heavier above the knee. However, results contradicting this conclusion have to be considered and understood. 



1. Introduction 

The energy spectrum of cosmic rays (CR) spans 
a very wide energy range, with particle fluxes 
steeply falling more than 30 orders of magnitude. 
Above the solar modulation region, the spectrum 
can be well described by a power law, which 
steepens around 3 lO^^eV, a feature called the 
"knee", discovered in 1958 ||]; it softens again at 
~ lO^^eV, the "ankle" . 

Explaining the knee feature would shed light on 
the CR origin and acceleration mechanism, de- 
pending on whether it is a signature of a change 
in the hadronic interactions at such energies or it 
reflects a feature of the cosmic ray spectrum, thus 
concerning mainly astrophysics. 
Several arguments involving energetics, composi- 
tion and secondary 7 ray production suggest that 
cosmic rays at least up to the knee region are con- 
fined in the Galaxy. 

The most popular theory is that of diffusive shock 
acceleration in Supernova remnants (SNR), that 
is particle acceleration by SNRs expanding super- 
sonically in the surrounding medium. 
Supernova explosions can easily account for the 
energy stored in galactic cosmic rays; the spec- 
trum emerging from the SNR is of the type E~^-^ 
up to a maximum energy near 10^* eV times the 
nuclear charge, after which it drops very rapidly. 



Folding the production spectrum with the effect 
of diffusion through the Galaxy, and taking a 
trapping time varying as E~^'^ (as found from 
the proportion of secondary to primary nuclei ar- 
riving to Earth), the resulting flux of CR in the 
Galaxy would be oc i?^^-^, in close agreement 
with expectations. The maximum achievable en- 
ergy is close to the knee one. 
A direct evidence that the nucleonic component 
of CR is indeed produced in SNRs could be ob- 
tained by the observation of 7 rays: the accel- 
erated cosmic rays can in fact interact with the 
local interstellar matter, in this way producing 7 
rays by either hadronic or leptonic production. 
Various experimental groups reported on TeV 7 
emission from supernova remnants, like SN1006, 
RJX1713. 7-3946, Cassiopea-A (see e.g. [g and 
references therein) . Unfortunately in all cases the 
emission can be attributed to electron progenitors 
and no positive evidence for hadroproduction of 
TeV 7's has been found yet. 
Further information can be obtained from the 
study of the distribution of CR arrival directions. 
A recent compilation of the anisotropy measure- 
ments can be found in |^; while the amplitude 
and phase of anisotropy data below ~ 2 lO^^eF 
are consistent and statistically accurate the 
experimental results at higher energies are not 
compatible with expectations. This could mean 
that the diffusion model cannot be simply extrap- 
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olated to these energies 

Various models have been put forward trying to 
identify the sites and mechanisms of injection of 
cosmic rays at higher energies, at and above the 
knee. If the bend in energy spectrum is related to 
the maximum achievable energy in the accelera- 
tor, then CR at higher energies could be powered 
by a reacceleration by interstellar turbulence ||], 
or they could be produced by Supernovae explod- 
ing in denser media (their stellar wind cavity) (Q 
and references therein). On the other hand, the 
knee could be attributed to propagation effects. 
In both cases, one would expect multiple bends 
due to the different elements bending at fixed 
rigidity; the composition would become heavier 
above the knee. An extra-galactic origin for CR 
above the knee has also been proposed Q, where 
the accelerator sites are found in Active Galactic 
Nuclei and the resulting composition is getting 
lighter above the knee. 

A completely different point of view assigns the 
knee to a new dramatic process of hadronic in- 
teraction which takes over around the knee en- 
ergy. However, even if it is true that in the AE 
of interest we have no direct information about 
the hadronic interaction cross section for the sec- 
ondary production relevant to the interpretation 
of measurements, no experimental data as far 
show a need for a different interaction mechanism. 
From the experimental point of view, what is 
most important in order to test the models is to 
measure the cosmic ray composition and energy 
spectrum near the energy limit of the shock mod- 
els; moreover, measurements of anisotropy and 
secondary to primary ratio at higher energy are 
of utmost importance. 

2. Direct measurements 

Direct measurements of the relative abun- 
dances of the cosmic ray nuclei and their distri- 
bution in energy are possible only at relatively 
low energy: they in fact require installation of 
instrumentation on balloons or space shuttles fly- 
ing outside the atmosphere at very high altitude. 
The most recent results still come from two bal- 
loon experiments, JACEE and RUNJOB, as sum- 
marised in |p|,p^. 



The proton and Helium spectra have been mea- 
sured by JACEE up to about 800 TeV |ll|; no 
break was found in the proton one, but above 
80-90 TeV the experiment does not have enough 
statistics to either assess or reject its presence 
111- 

JACEE claim for a flatter He spectrum as com- 
pared to the proton one is in agreement with 
SOKOL result but this is not confirmed 

by the results by RUNJOB The JACEE 

group reported values are 7p = (—2.80 ± 0.04) 
and 7ffe = (-2.68 ±0.06), while RUNJOB slopes 
are both ~ —2.80 with an uncertainty between 
10 and 20% for both protons and Helium nuclei. 
It should however be mentioned that the signif- 
icance of the difference between the slopes for p 
and He is at the level of only 2<t; on the other 
hand, the results from RUNJOB are based only 
on the 1995 and 1996 data (~ 30% protons and 
~ 13% He respect to JACEE data). The experi- 
mental results on the p and He slopes are of par- 
ticular importance as regards the models of non 
linear acceleration of cosmic rays, where the injec- 
tion rate is an increasing function of the primary 
particle rigidity |p^ . 

The single component spectra are shown in Fig.|l| 
and in Fig.^. The experimental results agree for 
what regards the iron spectrum, while RUNJOB 
gives a factor of 2 lower spectra for the C-N-0 
and Ne-Si groups. Data are all consistent with 
an increase of the mean logarithm of the average 
primary mass < In A > with energy, as shown in 
Fig.|. 

It is clear that more statistics is needed above 
100 TeV, and various new projects are under 
development. ACCESS ||l^ is estimated to be 
launched on the International Space Station in 
2006. Its primary goal will be the measure of 
energy spectrum and composition up to lO^^eV, 
thus testing Supernova shock acceleration mod- 
els; the charge range at high energy will be 1 < 
Z < 28. Three different detectors are being built 
for this purpose: a charge identification module, 
to measures the abundances of all individual ele- 
ments, a transition radiation detector, to identify 
and measure the energy of particles with Z > 2 up 
to ~ 100 TeV/nucleon, a calorimeter to measure 
the particle energies and to identify electrons. 



The ATIC project, in its initial design for 
long duration balloon flights, is devoted to study 
the energy spectrum of Galactic proton and he- 
lium up to lO^^eV, in order give information 
about the proton/helium ratio, the possible dif- 
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Figure 1. Differential energy spectra for proton 
and helium llOjl - 
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Figure 2. Differential energy spectra for the CNO, 
NeMgSi and Iron groups 
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Figure 3. Average primary mass vs primary en- 
ergy from direct measurements \1L 



ference in their spectral slopes, the existence of a 
bend in t he p roton spectrum. 
CREAM plans to explore spectrum and com- 
position up to ~ lO^^e^, exploiting ultra long 
duration balloon flights (~ 100 days). With 
an exposure of ~ 300 m?sr days, this instru- 
ment will collect ~ 500 proton and helium nuclei 
above lO^^'eT^, reaching ~ 30% statistical accu- 
racy above lO^^eV. 

The new Ionization Neutron Calorimeter INCA 
figj ] proposes to study the range 0.1-10 PeV us- 
ing the well known techniques of ionization and 
neutron monitor to measure energy and a silicon 
particle charge detector to determine the charge 
and coordinates of the primaries. 
Combining different detectors, the new projects 
will have a very powerful tool to overcome the 
individual technical limitations. 

3. The energy spectrum 

The experimental observables which are mea- 
sured in order to extract information about the 
energy spectrum are the charged components of 
showers as measured by ground based detectors 
with scintillator counters, muon and hadron de- 
tectors, or the Cerenkov light produced by shower 



particles as they propagate through the atmo- 
sphere. 

The interpretation of these ground level observa- 
tions in terms of primary particle characteristics 
is far from straightforward, being strongly depen- 
dent on models simulating the production and 
propagation of particles through the atmosphere. 
Models in turn depend on extrapolations applied 
to the data on high energy particle interactions 
studied at accelerators; the energy region of in- 
terest is in fact much higher than that studied at 
accelerators, the explored kinematic region is the 
forward one, the collisions among nuclei make the 
influence of nuclear effects not negligible. 
The most recent results of EAS experiments con- 
cerning the primary energy spectrum of cosmic 
rays are described in the following. 
The electron size spectra as measured by the 
EAS-TOP experiment [|0| are shown in Fig.|. 
The knee is clearly visible, and the size corre- 
sponding to the knee shifts towards lower values 
at increasing atmospheric depth, as expected for 
a knee at given primary energy ]2l| ] . 

The shower size at the knee decreases with in- 
creasing atmospheric depth, with an attenuation 
length Ak = (222 ±3) g cm~^, in very good agree- 
ment with that found for the shower absorption in 
atmosphere; the integral intensities /fc(> Ek) are 
constant within 20%. The knee in electron size 
is quite sharp, showing that the change in slope 
occurs in a limited range of {ANe/Ne < 25%). 

The muon size spectra have been measured in 
4 different zenith angle intervals, as one can see 
in Fig.||, where the change of slope is visible at 
all atmospheric depths despite the large statisti- 
cal fluctuations. The knee, around Nj^ ~ 10^'^^, 
is independent on the number of detected muons, 
being in fact visible at any core distance ||2^ . The 
integral fluxes in electron and muon size are com- 
patible at all atmospheric depths, as expected for 
a feature occurring at fixed primary energy, also 
confirming the consistency of the whole proce- 
dure. 

A further interesting result comes from the rela- 
tion between the electron and muon size slopes, 
which can be written as N^^ oc iV" with a ~ 0.75 
in all angular bins: no sudden change in the sec- 
ondary production when going through the knee 
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Figure 4. Differential electron size spectrum at 
different zenith angles, that is different atmo- 
spheric depths, as measured by EAS-TOP. The 
solid lines show the results of the fitting procedure 
Ji^; the fits represented by dashed lines were ob- 
tained requiring a constant integral flux above the 
knee. 

region is seen, thus showing that, at least from 
this point of view, no new hadronic effects is 
needed to explain the knee 

A simulation of the shower production and devel- 
opment in atmosphere using the CORSIKA code 
with the HDPM interaction model allows to 
find the relation between shower size and primary 
spectrum Ne{Eo,A) = a{Aeff )Ei^^^''"^ The ef- 
fective mass Aeff is calculated from the extrap- 
olation of the single nuclear spectra measured at 
low energies by direct measurements; above the 
knee, a rigidity dependent cutoff is used. 
The final result is shown in Fig.^; the agreement 
with direct measurements at low energies and 
with other air shower experimental results at the 
highest ones is quite good. The systematic uncer- 
tainties in the energy spectrum are due to the pri- 
mary composition and interaction model chosen 
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Figure 5. Differential muon size spectrum at 4 
different atmospheric depths as measured by EAS- 
TOP I? 



Figure 6. Primary energy spectrum as obtained 
by EAS-TOP compared with other experimental 
results /Ezl/. 



in the calculation: the maximum difference in the 
determination of Ne between different interaction 
models and HDPM is ~ 10%; the all-particle flux 
obtained with "heavy" or "light" limit composi- 
tions differs by ~ 10% from the above calculated 
one. 

The energy spectrum is determined by CASA- 
MIA using the muon and electron N* size 
measurements . The N* indicates in this case 
the sum of e+,e~,7 at the ground. The sizes 
combination F = logiQ{N* -I- ipNf^) was found to 
be log-linear in Eq and, what is most important, 
independent on the primary mass. 

In Fig.0 this relation is shown as found from a 
simulation of primary protons or iron; the model 
used was the QGSJET one. The systematic dif- 
ferences in energy assignment for different pri- 
mary mass A are less than 5%. The average abso- 
lute energy reconstruction errors go from 25% 
at lO^^eV to ~ 16% at > lO^^eV. 
The parameter ip, which defines the relative 
weight of muons and electrons in the showers, 
is strongly dependent on the model used for 
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Figure 7. F — logio{N* -\- ipNfj,) as a function 
of energy for simulated proton ( open circles ) and 
iron (black triangles) primaries |^/ . 



hadronic interactions, but the change in energy 
assignment due to this effect is claimed to be 
< 10%. The mass insensitivity allows to deter- 
mine the energy free of systematic effects (on the 
contrary, for example, if in some region the en- 
ergy spectrum changes, Ne vs E also changes). 
The energy spectrum thus derived is shown in 
Fig.||; the knee is located at the same primary 
energy for any atmospheric depths, as expected. 
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Figure 8. Primary energy spectrum as found by 
CASA-MIA M 



KASCADE [g^ measures all the three charged 
components of extensive air showers: electrons, 
muons and hadrons. 

The knee is clearly visible in all components; for 
electrons and muons, the higher statistics allows 
to study the size at different atmospheric depths, 
thus finding that the size at the knee decreases at 
increasing atmospheric depth |p7| . 
The primary energy spectrum can be extracted 
from the measured size spectra depending on the 
knowledge of the mass composition as obtained 
from the observables under investigation and on 
the relation between size and energy resulting 
from simulation. 

The energy spectrum shown in FigJfl was found 



by a combined minimisation to fit both the A^e 
and the A^^ truncated muon size spectra simul- 
taneously (the truncated muon number is that 
found by fitting the muon lateral distribution 
within a hmited range of 40-200 m) [|§. The 
evaluated size spectra are in fact the convolu- 
tion of the energy spectrum and a kernel function 
describing the probability of a given primary to 
produce a shower with a certain size and which 
includes the parametrisations of shower fluctua- 
tions for both proton and iron primaries accord- 
ing to Monte Carlo. The knee is found at about 
4 PeV. 
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Figure 9. Primary energy spectrum by KAS- 
CADE from electron and muon size spectra 



In FigjT^, the energy spectrum as derived from 
hadronic data is also shown. The knee is again 
very clear, and the expectations from pure beams 
of protons or iron primaries are shown for com- 
parison. Furthermore, KASCADE data show ev- 
idence of the knee also in the energy sum of 
hadrons in the calorimeter [p9| . 
From muon density measurements in the multi- 
wire proportional chambers below the central de- 
tector [Q, a subdivision of the data in "light" 
and "heavy" samples (according to the parame- 
ter Log N^/Log Ne as described in Sect. 4) shows 
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Figure 10. Primary energy spectrum as found us- 
ing KASCADE hadrons. Empty circles and trian- 
gles are the spectra from simulation of pure iron 
or proton primaries respectively \2!\]. 



that the knee is strongly dominated by the light 
component, within a 30% uncertainty due to 
Monte Carlo statistics. 

The broader lateral distribution of the Cerenkov 
light, due to the smaller absorption of photons in 
atmosphere, and the high photon number density, 
that means a better signal-to-noise ratio even for 
smaller arrays, are the main advantages in using 
Cerenkov detectors as compared to charge parti- 
cle counting arrays. The most recent results from 
apparata based on the detection of Cerenkov light 
from Extensive Air Showers come from BLANCA 
31 1 and DICE |^^, both operating at the same 
site and sharing some equipment with CAS A. 
The first one consists of 144 angle-integrating 
Cerenkov light detectors located in the CASA 
scintillator array, which provides the trigger and 
gives core position and shower direction. The 
Cerenkov lateral distribution function is mea- 
sured and fitted through the expression C(r) — 
Ci2oe~'^^ in the inner part of the distribution 
(r < 120 to). The intensity at a critical radial 
distance of 120 m, entirely determined by den- 
sity and scale height of the atmosphere, is propor- 
tional to the primary energy and the dependence 
on the primary mass is fully included in the slope 
s of the distribution, which is in fact a function 
of the depth of maximum development X^ax- 



DICE consists of 2 imaging telescopes of 2m di- 
ameter. What is measured is the Cerenkov light 
size N^, by summing the total amount of light at 
each phototube and the depth of maximum devel- 
opment of the shower X^ax by fitting the shape 
of the light image in each telescope, with a pro- 
cedure that is essentially geometrical and not de- 
pending on simulations, except for calculations to 
determine the angular distribution of light around 
the axis. The core position and shower direction 
are given by the CASA scintillator array. 
The primary energy is estimated through a fit in- 
cluding geometry, and X^ax and takes there- 
fore into account the dependence of the lateral 
distribution and intensity of the Cerenkov light, 
at fixed primary energy, on the primary mass. 
The resulting energy spectra from BLANCA and 



DICE are shown in FigjTl] (QGSJET model) 
and Figjl^ in comparison with other experimen- 
tal results. The knee feature is evident in the 
BLANCA energy spectrum at ~ 3 PeV] a 10% 
shift in the energy scale, which is however less 
than the instrumental uncertainty, is found by 
changing the interaction model chosen to inter- 
pret the data. The absolute calibration error re- 
sults in a ~ 18% systematic error on the energy 
assignment. According to DICE data, the knee is 
found around 3 PeV; the systematic uncertainty 
in the absolute flux is ~ 30%, due to the intrinsic 
error in the energy scale of ~ 15%. 

The Cerenkov lateral distribution has also been 
measured by HEGRA with its AIROBICC de- 
tectors iQ. Again, the scintillator array pro- 
vides the core position and arrival direction of 
the showers. The total primary energy can be 
reconstructed by the electromagnetic one, if one 
assumes a primary composition; however, a de- 
termination of the primary energy in a mass in- 
dependent way can be obtained following the ap- 
proach of Lindner , at the expenses of the en- 
ergy resolution, which is worse than that found 
by the mass-dependent method and of a stronger 
dependence of the result on the fluctuations of 
Xmax- The knee is found at about 3 PeV. 
In Fig.^, the all-particle primary energy spec- 
trum as obtained with the various experiments 
here described is shown. A summary of the up-to- 
date situation is given in Table |l|. The differ ences 
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Figure 11. Primary energy spectrum from 
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Figure 12. Primary energy spectrum from DICE 
events in coincidence with CA SA -MIA ( |?^/ ) . 

among the quoted knee energies are mainly due to 
the assumed composition, which in turn depends 
on the observables used as will be discussed be- 
low, but the existence of the knee is clearly estab- 
lished between 2 and 5 PeV. 



Figure 13. All-particle primary energy spectrum 
from the various experiments as described in the 
text. 

4. Composition 

The cosmic ray primary composition measure- 
ments around the knee are crucial for the under- 
standing of the mechanisms of acceleration and 
the source problem. 

The experimental observables we are dealing with 
are: a) the Cerenkov lateral distribution or the 
image of the Cerenkov light emitted by the 
shower in atmosphere; they allow to determine 
the depth of shower maximum Xjnaxj which is 
a logarithmically increasing function of the pri- 
mary energy. At fixed Eq, heavier primaries are 
expected to interact earlier, thus giving a smaller 
value of Xmax (higher in atmosphere), b) the 
muon and electron sizes of the showers. For 
a given primary energy, EAS induced by heavy 
primaries develop earlier in atmosphere and less 
energy is released in the electromagnetic com- 
ponent, thus producing a smaller Ne at ground 
level, as compared to proton showers; on the other 
hand, muons are produced more copiously in EAS 
by heavy primaries, because of the higher number 



Table 1 

Slopes of the energy spectrum below (ji) and above (^2) the knee, 
seen (PeV). 



Ek is the energy at which the knee is 
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2.76 ±0.03 


-3.19 ±0.06 


2.7-4.9 


2.70 ±0.05 


-3.10 ±0.07 


4.0-5.0 


2.66 ±0.12 
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smooth 
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-3.02 ±0.05 


-4.7 


2.60 ±0.04 


-3.00 ±0.05 


smooth 


2.60 ±0.02 
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~ -2.7 


~ -3.0 


~ 3.0 


2.67 ±0.03 
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of low energy pions. 

The EAS-TOP group studied the composition by 
analysing the behaviour of iV^j as measured in ver- 
tical direction in narrow bins of N^, correspond- 
ing to ANe/Ne = 12%. The result is shown in 
Fig.|l^: data are compared with the results of a 
full simulation including the detector response, 
where the 1 TeV composition with equal slopes 
for all components was used, in this way assign- 
ing constant composition with energy. 

The EAS-TOP data clearly suggest a growth of 
the mean A with energy, that is a heavier compo- 
sition above the knee. A change of ALog{Ne) — 
0.5 resuhs in a AA/A ~ 0.4 |§. 
The "K Nearest Neighbour" test was used by 
CASA-MIA to study the composition Us- 
ing the electron and muon densities at different 
distances from the core and the slope of the elec- 
tron lateral distribution, samples of event for each 
different primary mass are generated by Monte 
Carlo. An experimental event is assigned to the 
"light primary" or "heavy primary" class by look- 
ing at the K nearest neighbours (KNN) in the 
plane of the used variables: the event will belong 
to the light primary group if e.g. more than 50% of 
its KNN are light primaries. Due to fluctuations, 
which tend to superimpose classes, only broad 
classes of "p-like" and "Fe-like" events can be 
used. The proton resemblance, defined as the av- 
erage fraction of K nearest neighbours which are 
protons, is shown in Fig.|l^ for K=5, normalised 
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Figure 14. < A^^ > vs N,, from EAS-TOP experi- 
mental data (full squares ) as compared to a simu- 
lation with mixed composition and all components 
with the same slope (downward triangles), pure 
proton (upward triangles) and pure iron (stars) 
primaries / [^/ . 

such that a pure proton composition would lay 
along the top of the plot and a pure iron along 
the bottom border. 

The trend towards a heavier composition above 
the knee is evident; a change in the hadronic in- 




Figure 15. Normalised proton resemblance plot 
from CASA-MIA. The open squares give the es- 
timated result if the composition is taken from 
JACEE direct measurements 

teraction model used in the simulation does not 
change the result. Classifying the events accord- 
ing to their probability of being light or heavy 
primaries and plotting the energy spectrum for 
the two classes separately, CASA-MIA data sug- 
gest that the knee be due to the hght mass group; 
the spectra are consistent with the idea of cutoffs 
proportional to the particle rigidities. 
The composition problem has been attacked by 
KASCADE people in a variety of ways, using dif- 
ferent observables and analysis methods. 
The most sensitive dependence on primary mass 
was identified in the ratio log N^^ /log N^, which 
is found to be Gaussian distributed at fixed A 
. The experimental ratio is fitted by a super- 
position of simulated distributions (one for each 
primary mass group), directly obtaining the frac- 
tion of each mass group, as shown in Fig.|l^ for 
two energy bins. The composition is dominated 
by the light component up to about 4 PeV, get- 
ting heavier above the knee; the analysis also 
proves that the composition cannot be described 
by a single component. The < In A> so obtained 
is shown in FigJT^. 

A number of hadronic observables has also been 
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Figure 16. log Nj/" /log Ne for two of the consid- 
ered energy bins: 6.2 < log{E / GeV) < 6.3 and 
6.7 < log{E/GeV) < 6.8 

used, such as the lateral hadron distribution, the 
hadron energy spectrum, the maximum hadron 
energy etc. in order to investigate the composi- 
tion. As one can see in the world survey given 
in Fig.|l^, the hadronic data alone give a heavier 
composition as compared to other data. 
An interesting approach was used in ||4^ , where a 
multivariate analysis using all the measured com- 
ponents of the EAS is performed. The result 
shows a tendency to a lighter composition ap- 
proaching the knee, followed by an increase in 
the average mass above it. 

A comparison among the results by KASCADE 
shows that the absolute scale strongly depends 
on the observables which are used. It is clear 
that the balance of energy among the different 
components of EAS in the simulation does not 
reproduce the real situation; tests of the high en- 
ergy interaction models are been performed p4| . 
In experiments like BLANCA, Spase- VULCAN 
H, CACTI [|6), Hegra-AIROBICC, X^ax is 
measured from the slope of the Cerenkov lateral 
distribution, which is an almost linear function of 
the depth of shower maximum. This function is 
rather independent on the models chosen for the 



description of hadronic interactions, while any in- 
terpretation of the experimental results in terms 
of primary composition is not. 
In the case of DICE, the imaging technique allows 
to measure X^ax in a rather direct way, by fitting 
the shape of the shower Cerenkov image in each 
of the 2 telescopes, knowing the arrival direction 
and the core position of the shower. 
A survey of the results is shown in Fig.|7[ up to 
the Fly's Eye energies (where air fluorescence is 
measured); the "direct" point shows the X^ax 
that would be expected on the basis of balloon 
direct measurements ETT. 
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Figure 17. Mean height of shower maximum vs 
energy as measured by various devices. HEGRA 
data from Fly's Eye data from The 

lines show the expectations from a pure proton 
or iron composition using CORSIKA + QGSJET 



composition becoming progressively lighter with 
increasing energy. 

DICE and CASA-MIA groups studied the com- 
position problem also by means of a combination 
of measured parameters [Q. Two estimates of 
the mass have been derived, one using the X^ax 
as determined by DICE and the other with A^^ 
and Ne by CASA-MIA, for each detected shower. 
The combined use of different measurements al- 
lows first of all a study of the systematics biasing 
the composition results; moreover, the require- 
ment of consistency among various measurements 
allows to limit the range of the parameters used 
in the models. This analysis suggests a primary 
composition becoming lighter at and above the 
knee, not excluding however a constant composi- 
tion around the knee energy. 
A survey of the previously described results on 
primary composition in terms oi < In A > \s 



shown in Fig, 18 
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BLANCA data suggest a composition getting 
lighter near the knee and turning to a heavier one 
after the knee energy. Data from DICE require a 



Figure 18. Mean logarithmic mass vs primary en- 
ergy. The box represents the region of direct mea- 
surements ^1/. MSU data from Chacaltaya 
data from p3J. 



5. Conclusion 

The cosmic ray energy spectrum and compo- 
sition are studied with a variety of experimental 
techniques detecting different air shower compo- 
nents in the energy region above 1 TeV. 
Below the knee, no new data are available and 
the conclusions reached by JACEE and RUNJOB 
experiments still hold; however new projects, 
planned to fly on balloons or on the Space Sta- 
tion, are in progress. They will surely extend the 
explorable energy region and the available statis- 
tics on single nuclei. 

The energy spectrum has been studied in de- 
tail both by charged particles and Cerenkov light 
ground arrays and some firm conclusions were 
reached: all data agree on the existence of the 
knee in the primary energy spectrum of cosmic 
rays at an energy ~ 3 — 4 PeV. The bend has 
been seen in all shower components, thus support- 
ing an astrophysical interpretation of the knee as 
opposite to that of a change in the hadronic inter- 
action picture at these energies. All results agree 
to attribute the knee to the medium-light mass 
primaries. 

Progress has been made as regards the mass com- 
position. Almost all the ground array results 
show an increase in the primary mean logarithmic 
mass above the knee, even if the absolute scale 
can be quite different. There are however contra- 
dicting results coming from experiments relying 
on the Cerenkov light detection from Air Show- 
ers; the differences in the measure of Xmax are 
however quite big. 

It is very important to study in detail the system- 
atics which could bias the results; the combined 
use of different observables and the comparison 
among various data sets can help in this task. The 
different sensitivity to composition of the various 
used observables, the methods employed to deter- 
mine the primary energy and the problems found 
in the models used in the simulations could ex- 
plain the spread in the results which is apparent 
in Fig.|l|. 
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